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ABSTRACT 


The lens of the eye has no blood vessels to interfere with vision. The lens is far 
too large for diffusion to provide food and clear wastes. 

Experimental, theoretical and computational work has shown that the lens 
supports its own microcirculation. It is an osmotic pump that implements what 
physiologists have long believed “Convection provides what diffusion cannot.” 

We introduce a general (non-electro-neutral) model that describes the steady- 
state relationships among ion fluxes, water flow and electric field inside cells, 
and in the narrow extracellular spaces within the lens. 

Using asymptotic analysis, we derive a simplified model exploiting the numerical 
values of physiological parameters. The model reduces to first generation 
‘circuit’ models and shows the basis of computer simulations too large to easily 
understand. The full model helps resolve paradoxes that have perplexed 
molecular biologists: crucial physiological properties do not depend as 
expected on the permeability of the lens interior (to water flow). 
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Theme for Recent Lectures at Kansas U 

Bob Eisenberg 


Examples of Precise Biophysics 
at different scales, that yield 
Physical Understanding of Biological Systems 

and How They Work 


rV>e cV^ ed .^ e v' la ' : ' 0< ' 


The Lens of the Eye: 

Bidomain Model of an Osmotic Pump 



1 cm 
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The Lens is a Device to gather and focus light. 
It is designed by evolution for a purpose: 

To gather and refract light 








The Lens gathers and refracts light. 
No blood vessels!!! 
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The Lens is a 

Preparation 


Chosen to Study 
Convection and Circulation 
By an Osmotic Pump 


Without blood vessels! 



Surface venules 


urface arteriole network 


Penetrating venules 


Columnar boundary 


Subsurface 

microvessels 


Penetrating arteriole 


Circulatory System is Incredibly Dense. 
Convection is needed. 

Convection delivers what 
Diffusion cannot! 

Diffusion is too slow to support life!!! 





1 cm 


Central Question: 

How does the lens support circulation 
without blood vessels? 







Central Idea of Biological Experiments 

Preparation 

We imagine that other biological systems 
use the same mechanism as the lens, 

most importantly 

the tiny extracellular spaces in the brain 


Central Question: 

How does the lens support circulation 
without blood vessels? 
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Central Question: 


How does the lens support circulation without blood vessels? 

ANSWER: 

The Lens is an 

Osmotic Pump 








The Lens 


An Osmotic Pump 

Can we prove that the lens is an 
Osmotic Pump? 


Are Osmotic Pumps a 
GENERAL BIOLOGICAL MECHANISM? 

This is a seminar on Quantitative Biology as well as math 


Crystalline lens 


Equator 

Anterior 

surface 

Posterior 

surface 


Epithelial ... 

Lens 

nucleus 

"a 


Capsule 

Cortex 




*■ 


1 r. m 
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Start with the Electrical Structure 

because 

Electrical Methods are Most Precise 


Structural Analysis of Electrical Properties 


1. Qualitative Structure = Anatomy 

2. Quantitative Structure Stereology 

3. Electrical Model (from Poisson) 

4. Impedance Measurements (DC to 10 4 Hz) 

5. Parameter Estimation 

6. Checking Parameter Invariance 

7. Prediction of Biological Function 
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Structural Analysis of Electrical Properties 

uses 

MATHEMATICS 


1. Qualitative Structure = Anatomy 

2. Quantitative Structure Stereology 

3. Electrical Model (from Poisson) 

4. Impedance Measurements 

5. Parameter Estimation 

6. Checking Parameter Invariance 

7. Prediction of Biological Function 


INVERSE PROBABILITY 
FIELD THEORY 
CIRCUIT THEORY 
INVERSE PROBLEM 
INVERSE PROBLEM 
FIELD THEORY 
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Electrical Methods are Most Precise 
Spherical Cell 



Electrical Models often using 

Singular Perturbation Theory 

exploiting membrane boundary condition 

1 .J Gen Physiol, 1970 . 55(6): p. 736-57. 

2. Prog. Biophys. Mol. Biol., 1970 . 20: p. 1-65 

3. SIAM J. Appl. Math., 1971 . 21(2): p. 339-354. 

4. Biophys J, 1972 . 12(4): p. 384-403. 

5. Anna Rev Biophys Bioeng, 1973 . 2: p. 65-79. 

6. J. Mathematical Biology, 1975 . 2: p. 277-300. 

7 . J. Mathematical Biology, 1975 . 2(4): p. 301-316. 

8. SIAM J. Appl. Math., 1976 . 30(2): p. 222-239. 

9. Ann N YAcadSci, 1977 . 303: p. 342-54. 

10. Biophys J, 1977 . 17(1): p. 57-93. 

11. Biophys J, 1978 . 23(2): p. 277-84. 

12. CRC Crit Rev Bioeng, 1980 . 4(3): p. 203-32. 

13. Biophys J, 1983 . 44(2): p. 225-48. 

14. Biophys J, 1983 . 42(1): p. 55-9. 

15. “Membranes, Channels, and Noise 1984”, p. 49-116. 

16. Biophys J, 1985 . 48(2): p. 253-67. 

17. Kevorkian, J. and J.D. Cole, 1996 

Multiple Scale & Singular Perturbation Methods. 
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Electrical Models of Syncytia 
Can Be Derived Using Definitions of Vector 
Operators, Div and Grad 



Eisenberg, Barcilon, and Mathias. Biophys. J. 25: 151-180 (1979). 





















The Lens: An Osmotic Pump 

Electrical Methods are Most Precise 
Impedance Measurements of Spherical Syncytia 


CURRENT PASSING 



Mathias, Rae, and Eisenberg. Biophys. J. 34: 61-85 (1981). 
























































Electrical Measurements are Precise 
Spherical Syncytia 
TYPICAL DATA not selected for fit 



Figure 9 The magnitude, phase angle, and theoretical fit to impedance data of lens 2-8 taken from one 
electrode location. The electrode was at r, = 0.93 mm in a lens of radius a = 1.33 mm. The circuit 
parameters determined from the outermost location of the electrode are quite close to the average 
parameters determined from many locations (Table IV). 


Rae, Mathias, and Eisenberg. Exp. Eye Res. 35: 471-490 (1982). 
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How is this precision possible? 
Biology has Invariants ! 


Cell Interior is a (nearly) pure resistor 
Membrane is a (nearly) a pure capacitor 
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Biology has Invariants! 


Cell Interior 

is a (nearly) pure resistor 


. _ v 

1 ~ R 
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How is this precision possible? 
Biology has Invariants! 

Interior of Cells is (nearly) a pure resistor 

even a muscle fiber 
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PHASE ANGLE 


Muscle Interior is a Pure Resistor 

Biology has Invariants 

TABLE I 

RESULTS FROM 14 MUSCLE FIBERS 


0max 


degrees 

-1.48 

1.33 

0.74 

-2.45 

-2.42 

- 2.22 

-1.58 

-1.60 

0.59 

0.68 

-1.79 

-0.62 

- 1.86 

-2.25 

Mean —1.07 

SD 1.34 





o o 


O o 


OOOOOOOOOOOOOOOOOO 


1 10 100 1,000 10,000 

LOG FREQUENCY (Hz) 22 


Mobley, Leung, and Eisenberg, J Gen Physiol, 1975. 65: p. 97-113. 






Interior of Cells is (nearly) a pure resistor 

even a muscle fiber 
How is this possible? 

Electric Current Flows AROUND Obstacles 

Electric Current Flows through Low 

Resistance Path 

Property of PARALLEL RESISTORS 
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Biology has Invariants! 


Cell Membrane 
is a (nearly) pure capacitor 

dV 

i = C — 
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Glycoprotein: protein with 
carbohydrate attached 


, Glycolipid: lipid with 
/ carbohydrate 
attached 


Peripheral membrane 
protein 


Integral membrane ^ 
protein 


Cholesterol 


Channel protein 


Phospholipid 

bilayer 


Denis Haydon Life Work 

Membrane is a (nearly) a pure capacitor 

1. J Theor Biol, 1965. 9(3): p. 433-43. 

2. J Theor Biol, 1965. 9(3): p. 422-32. 

3. J Theor Biol, 1965. 9(2): p. 278-96. 

4. J Gen Physiol, 1965. 48: p. SUPPL:59-63. 

5. J Mol Biol, 1968. 32(1): p. 149-50. 

6. Br Med Bull, 1968. 24(2): p. 124-6. 

7. J Theor Biol, 1968. 18(3): p. 371-9. 

8. J Am Oil Chem Soc, 1968. 45(4): p. 230-40. 

9. J Membr Biol, 1971. 5(3): p. 277-96. 

10. Biochim Biophys Acta, 1979. 557(1): p. 259-63. 

11. J Physiol, 1979. 287: p. 38P-39P. 

12. J Physiol, 1979. 287: p. 2P. 

13. Biophys J, 1980. 30(1): p. 129-36. 

14. J Physiol, 1980. 309: p. 229-45. 












Estimated * Membrane Capacitance 

of Lens 
does NOT vary 

with ion concentrations 


Expanded 

Scale 



BATH POTASSIUM CONC (MM) 


Fig. 3. Plot of the capacitance ( ' m of t he inner membranes vs. the potassium concentration in the bath. 
The capacitance was determined by fitting the syncytial theory of Eisenberg et al. (1979). to the lens 
impedance data. Values are the mean + s.n. with the number of lenses measured at each concentration 
given in parentheses. The mean capacitance determined from all measurements is shown (-). 


*Note: this estimation involves structure, qualitative and 
quantitative, impedance measurement, field theory, and 
solving inverse problem. Errors in any of these will 
produce estimates that vary with concentration. 
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The Lens 


An Osmotic Pump 


Can we prove that the lens is an Osmotic Pump? 
Life work of Richard (Rick) Mathias 



State University of New York SUNY Stony Brook 


Too many papers to list! Do Google Search 
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3 iuo/s Z * 



Mathias, Kistler, and Donaldson, 

J Membr Biol, 2007. 216(1): p. 1-16. 


28 


Na + Leak Conductance 







Computer Model 

of Lens 



B 




ion channels 

Na«7K+ 

ATPase 

pumps 


K+ channels AQpl 

Gap junctions AQpO 


Vaghefi, E., N. Liu, and PJ. Donaldson, Biomed Eng Online, 2013. 12: p.85 
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Direct Measurement of Pressure in the Lens 




Figure 2. The intracellular pressure measuring system. (A) A sketch of the effect of intracellular pressure on the electrode-intracellular 
solution interface. An elevated intracellular pressure (e.g., 100 mmHg) over that in the bathing solution causes the interface to move up 
the shank of the electrode, thus filling the narrow tip with relatively high resistance cytoplasm. (B) When the same 100-mmHg pressure 
is applied to the port of the microelectrode, the interface moves back to the tip as shown, and the electrode resistance is restored to its 
relatively lower value recorded in the bathing solution. (C) The manometer used to adjust the hydrostatic pressure at the electrode port. 
The crank drives a piston to create the pressure, which is connected to the electrode port through the plastic tubing. When the electrode 
resistance is restored to its value in the bathing solution, such that an increase in applied pressure has no effect on resistance, whereas a 
small reduction in pressure causes a small increase in resistance, we assume the applied pressure equals the intracellular pressure. The 
value of applied pressure is dien read from the column of mercury. 

1 A Ant Dnlo , . C 


Gao, Sun, Moore, White, Brink, and Mathias, 
J Gen Physiol , 201 7. 137: 507-20. 



Gap Junction Channels 


r/a 























Theory 



Gao, Sun, Moore, White, Brink, Mathias 
J Gen Physiol, 2011 
137:507-20. 


Direct Measurement of 
Pressure in the Lens 


Experiment 


















Normalized p, 


Pressure with Reduced Transport in the Lens 



Normal Transport 


Reduced Transport 


Gao, Sun, Moore, White, Brink, Mathias 
J Gen Physiol, 2011 
137: 507-20. 
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Changing Cell to Cell Coupling 
by Molecular Genetics 




Increasing the Number 
of Junction Proteins 


Reducing the Number 
of Junction Proteins 


Gao, Sun, Moore, White, Brink, Mathias 
J Gen Physiol, 2011 
137:507-20. 
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p. 509 of Gao, Sun, Moore, White, Brink, and Mathias, 
J. General Physiology, 2011,137: 507-520 


“A comprehensive model of lens transport, 
including osmotic 

and hydrostatic pressure gradients in both 

the intracellular 

and extracellular compartments of the lens, 
is beyond the scope of this paper.” 



Mathias, Rae, and Baldo (1997) Physiological Reviews 21-50. 



Our Contributions 


are to Derive Consistent Field Equations 

and Solve Them 

Multiple Fields are Hard to Deal with 
Consistently without Variational Methods 

(in engineering models or computations) 

A great deal of confusion can result 
from nonmathematicai reasoning 
Including paradoxes 
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Osmosis through a Semi-permeable Membrane 

a Consistent Approach to Interactions 

arXiv:1806.00646 



Shixin Xu Zilong Song Huaxiong Huang 

A Bidomain Model for Lens Microcirculation 

Biophysical Journal (2019) 116: p. 1171-1184 





Yi Zhu Shixin Xu Huaxiong Huang 
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Inconsistent Models 

produce confusion 
even when correct 



replaces 

inconsistent Models Models 
with Consistent 

Partial Differential Equations 
and boundary conditions 
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Scientists and Poets can Reach 

but 

Engineers must Grasp 

and not just reach 
if 

Devices are to Work 

Uncalibrated Devices do not Work! 


Poets 

hope we will never learn the difference between dreams and realities 
“Ah, ... a man's reach should exceed his grasp, 

Or what's a heaven for?” 

Robert Browning 

"Andrea del Sarto", line 98 






Mathematics Creates our Standard of Living* 

Partial Differential Equations 

are needed to describe 
Flow 

driven by multiple fields 


*e.g., Electricity, Computers, Fluid Dynamics, Optics, Structural Mechanics, 
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We begin at the beginning: 


Ionic Solutions are Complex Fluids 

in which 

'everything interacts with everything else’ 

and 

Flows are driven by convection, and diffusion. 


Shixin Xu , Eisenberg, Song, and Huang 

( 2018 ) 

arXiv:1806.00646, 35 pages. 
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Bulk Solutions and Membranes 
are described by 
Consistent Analysis of 
Coupled Water and Ionic Diffusion and Flow 

The Energy Variational Principle and Sharp Boundary Methods 

have been applied to a variety of membrane models, 
allowing density of solutions to be a function of concentrations, 
as is seen every day in chemistry laboratories 

apparently for the first time 

Unfortunately, so far, analysis has only been performed for ideal solutions. 

Challenge 

No one seems to know 

how to formulate coupled water and ion flow and diffusion, 
in the realistic nonideal case found in biology and electrochemistry 

NOT equilibrium, not ideal, 
must deal with saturation and finite size ions 


Shixin Xu , Eisenberg, Song, and Huang (2018) 

arXiv:1806.00646, 35 pages. 
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Challenge 

No one seems to know 

how to formulate coupled water and ion flow and 

diffusion, 

in the realistic nonideal case found in biology and 

electrochemistry 

NOT equilibrium, not ideal, 
must deal with saturation and finite size ions. 


Suggestion: 

Combine EnVarA with Poisson Fermi approach 

of Jinn Liang Liu, 

more than anyone else 
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Starling’s Law has been derived: water flow is proportional to hydrostatic and osmotic 
difference from the Energy Variational Principle and a Sharp Boundary Methods 

Total Energy Functional . 

E tot = E kin + E int + E r 

~ E k in + Ees + Ei on + £V 

= Z £ + ( e kin + e es + ef on )dx + jT(e r + Yo)dS 3 

+ ((%[«’ + *>)* 

where c 0 is a characteristic ion density, C m is membrane capacitance, y () is the membrane surface 
tension, hi the following, square brackets always denote the jumps across the interface. 

£± fn± f ± dx = f n+ f + dx + / f~dx is the sum of integration in left and right compartments 


Dissipation IS NOT unique 
This is a PHYSICAL Model 
e.g. it must it time and 
concentration dependence and 
much more 
BUT 

Mathematics must be consistent 
with the 

TOTAL ENERGY FUNCTIONAL 


weighted by the corresponding value of /. 

The dissination functional is defined as 

A= Xf n+ 2 ^\ D ^ 2dx+ Z j n / ±lv ■ u±{2dx i+Z D t c i\ v p-r\ 2dx 


+ 


|G 1 ([/I £ ])d5 + ! 


(4) 


G 2 (Q P )dS 


109 


fli = 


Where D t} = (Vu + ( Vu) r )/2 is rate of strain, T] ± and A ± are the two Lame constants 
pp(cp, 0 ± ,p ± ) and D~ are the chemical potential and diffusion coefficient of i th ion. The first three 
tenns in dissipation functional are the dissipation induced by fluid friction, volume change and ion 
diffusion in the bulk region. The last two terms are the dissipation induced by irreversible osmosis .on 
the membrane. The forms of G X (x) > 0 and G 2 (x) > 0, for any x, will be discussed later. 



NOT unique 

This is a PHYSICAL Model 
e.g. it does not include 
Saturation produced by 
Finite size of ions 






Intracellular and Extracellular Velocities (1) 


The intracellular velocity depends on the gradients 
of hydrostatic pressure and osmotic pressure (41.51.62), and the extracel¬ 
lular velocity is determined by the gradients of hydrostatic pressure and 
electric potential (41,67), 

(2a) 


(2b) 


where 0/ is the electric potential in the / space, r c is the tortuosity of extra¬ 
cellular region, /j, is the viscosity of water, k e is introduced to describe the 
effect of electro-osmotic flow, and k / is the permeability of the intracellular 
region (/ = in) and the extracellular region (/ = ex), respectively. 


and 


Mex 


K ex d 


cl 


l . Pex k e T c (fiex 

/i dr dr 


Min 


K in ( d d 

— T Pin - yJ B T-o in 

fi \dr dr 








Flows and Coupling 


Then, we obtain the following system for intra- and extracellular veloc¬ 
ities in domain Q = [0, /?J: 


Coupling: 


2 / 0 ^^ex^ex) 

V CIV 

— —MyL m (P ex — P in + J m kgT (Oin — O ex )) 


(la) 


and 


Conservation: 


1 d 

^di- 


(/ {-AAex^ex 


+ AA in Ujn ) ) — 0 , 


(lb) 


where u t and P/ are the velocity and pressure in the intracellular and extra¬ 
cellular space, respectively, and 0 { is the osmotic pressure with definition 


Osmolarities: 


O 


ex 


= y ^c o- = y ^c 


in "b | 


•in 


in 


where C\ is the concentration of ith species ion in / space and A in /V in is the 
density of the permanent negatively charged protein. In this work, we as¬ 
sume the permanent negative charged protein is uniformly distributed 
within intracellular space with a valence of z. Here. .VI/ is the ratio of intra¬ 
cellular area (/ = in) and extracellular area (/ = ex), jM v is the membrane 
area per volume unit. y m is the intracellular membrane reflectance, L m is 
intracellular membrane hydraulic permeability. k B is the Boltzmann con¬ 
stant. and T is temperature. 
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Ion Circulation 

With similar assumptions, the conservation of ion concentration yields the 
following ion flux system: 

~2 J r = M,/ m (4a) 


+ = 0. (4b) 


The ion flux in the intracellular region J‘ in and ion flux in the extracellular 
region J l are defined as 
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cl . ■ 

C U e v - D‘ T c —C ~ D‘ T c 

ex ex ex c ex ex c K 




(5a) 


Membrane Flux of Ions 

j‘ m = ^ -</>«-£')’ (6a) 

fs = ^{<t>in-<t>ex-E‘), (6b) 

where E‘ = k^T/ e-\og{C‘ ex / C' m ) is the Nernst potential (an expression of 
the difference of chemical potential) of the ith species ion. 
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Electroneutral Expansion 

from PNP = Poisson Nernst Planck, 

named by Eisenberg & Chen, 1993, Biophysical Journal 64: A22. 

NOT electroneutral 


Because of the capacitance of the cell membrane, assumptions of exact 
charge neutrality can easily lead to paradoxes because they oversimplify 
Maxwell’s equations by leaving out altogether the essential role of charge. 
We use the analysis of (70) and thus introduce a linear correction term re¬ 
placing the charge neutrality condition (41.5 1) without introducing signif¬ 
icant error (see also (71 )), 


(1 - ij) ( ^>'<4, = —M r C m (4> in - 4> ex ) (10a) 


and 


A; 




= M v CMin ~ 0J, dOb) 

; * in / 


where i] is the porosity of intracellular region and C m is capacitance per unit 
area. 


Song, Cao, and Huang. Physical Review E, 2018. 97(1): p. 012411 
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Simplified Model 

by well-defined 
asymptotic expansion 


With ordering 

c 0 (e) 5 {^o 5 ^3} c 0(e _ ) 5 

{d 2 ,d^d 5 l d 6 ,d 7 }(Zo(e-). (is) 


d 9 = Df 1 — Df and <5 10 = Df 1 — Df cl , I = in, ex, it yields 

d 9 = 0(e 2 ), 5 I0 = O(e). (17) 


<*0 

<*1 

S 2 

<*3 

s 4 

$5 


M 


ex 


M in § 


k e T c k B T 

eRii* x 

fiRn 


in 


KmlrnksTO*' 

P* 

7 m k B TO*' 

M in ii* n 

RM v L m y m k B TO* 


u 


in 


Ls'Ysk B TO* 

M v C m k B T 


<*7 

<*8 


M v C m k B T 
e 2 C*( 1 — rj ) 
•Mex^ex 
MinD1 


in 


69 = 


<*10 = 
<*11 = 


D, c/ - Df 
Df ' 
D Cl _ jjNa 

w~~ 

r ci,o 

^in 

r cijo' 

^ex 


l G {in, ex} 
l G {in, ex) 















. (mmHg) 


Comparison with Experimental Results 
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A 500 


Dependence on ‘Connexin Permeability’ K in 



FIGURE 3 Comparison between different K in . (A) Comparison between simulation results and experimental results. The experimental data of dog, rabbit, 
and rat come from (47). Mice and Cx46 KI mice come from (28). According to (28), the Cx46 KI mice lens has twice the number density of lens gap-junction 
channels compared to mice. The parameter = 4.6830 x 10 20 /m 2 , and radius is written in dimensionless units for different species. (B) Space distribution 
of intracellular osmotic pressure (O in ). (C) Space distribution of intracellular electric potential (</>„,). (D) Space distribution of intracellular water velocity 
(»,,,). To see this figure in color, go online. 
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Non-Electroneutral Model vs. Simplified, and Mathias Models 
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Radius (mm) 



—Non-EN model 
-- Simplified model 





Mathias model 
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Resolution of Troubling Paradoxes 

in our Models 

Perturbation Expansion shows that 

Many intracellular and extracellular quantities 

are INSENSITIVE 
to increases in permeability 
EVEN BY A FACTOR of 20* 

Because many quantities 

depend only on the extracellular potential, 

and the 

Extracellular Potential is Insensitive 

to 

Permeability 
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Extracellular Electric Potential (mV) Extracellular Pressure (mmHg) 


Effect of Permeability K in 



Radius(r/R) 


Radius(r/R) 




Figure D.5: Comparison between different k ;m . 



































What next? 


“Brain Extracellular Space: The Final Frontier” 


Nicholson and Hrabetova, 
Biophysical Journal, 2017. 113: p. 2133-2142. 



Stirred by an Osmotic Pump 
like the lens, 
at work and in sleep ?? 
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Any Questions 



